CRPS-93-31 



c ^-. Comment on "Magnetic Susceptibility of the 

§5 : two-dimensional Hubbard model " 

CI,! Liang Chen and A. -M.S. Tremblay 

jjri ■ Departement de physique et 

Centre de recherche en physique du solide, 
Universite de Sherbrooke 
! Sherbrooke, Quebec, Canada, J1K 2R1 

O 
O 

On 
O 

m 

On 

: ABSTRACT 

a 

The observed magnetic spin susceptibility of high-temperature superconductors such 

O 

as La2- x Sr a; Cu04 increases when x increases from zero, i.e. as one dopes away from 
half-filling. Recent Monte Carlo simulations of A. Moreo (Phys. Rev. B 48, 3380 
(1993)) suggest that this behavior can be reproduced by the two-dimensional Hubbard 
model only at large coupling, namely, U/t of order 10. Using longer runs, our Monte 
Carlo simulations show that the same behavior as for U/t = 10 is obtained even in 
the intermediate coupling regime (U/t = 4), as long as the temperature is low enough 
(T = t/Q) that strong antiferromagnetic correlations are building up at half-filling. 
These results are consistent with the fact that in two-dimensions, the GRPA should 
fail in the parameter range where it predicts a magnetic phase transition. 
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One of the main puzzles in establishing the applicability of the one-band 
Hubbard model to high-temperature superconductors is the behavior of the uniform 
magnetic spin susceptibility as obtained, say by Knight-shift measurements. 1 While 
the maximum as a function of temperature is easy to understand qualitatively from 
the behavior of the Hubbard model at half-filling in the strong-coupling (Heisenberg) 
regime, the behavior away from half-filling and in the intermediate- or weak-coupling 
regime is theoretically much more uncertain. The limiting low-temperature value of 
the magnetic susceptibility experimentally increases as a function of doping, a behavior 
which at first glance can be reproduced within perturbative approaches only in models 
which include second- neighbor hopping 2 . With only nearest-neighbor hopping, non- 
perturbative effects, such as pseudo-gaps 3 are usually invoked. 

To numerically check under what conditions the one-band nearest-neighbor 
Hubbard model reproduces the experimentally observed behavior 1 , Moreo 4 has studied 
the magnetic spin susceptibility using determinantal quantum Monte Carlo techniques 5 ' 6 
She has shown that on a 4 x 4 lattice at inverse temperature f3 — 4/t the susceptibility 
has a maximum at half-filling and decreases with doping in the intermediate-coupling 
regime U/t = 4, while in the strong-coupling regime U/t = 10, the susceptibility 
increases with hole doping near half filling, resembling the experimental results 1 ob- 
tained for L^-z Sr^CuC^. She suggests that experimental results can be reproduced 
only when the coupling is quite large, namely U/t — 10 or more. 

We have used the same quantum Monte Carlo technique on a 4 x 4 lattice 
for U/t — 4 but at lower temperature, namely ksT = t/6, or inverse temperature 



[3 = 6/t. As shown in Fig. 1, we find that at this lower temperature, the susceptibility 
has qualitatively the same behavior as that observed for larger coupling (U/t = 10) 
and larger temperature {(5 — 4/t ) by Moreo. We measured the structure factor instead 
of the uniform magnetic susceptibility but these quantities are easily related using the 
fluctuation-dissipation theorem. To obtain the results of Fig. 1, one must perform 
much longer runs than is usually done. Typically, to obtain one point on the curve 
we have used At =1/15 and 100,000 warm-up sweeps of the whole space-time lattice 
before the measurements were taken. Measuring the static structure factor after every 
update of the lattice in space, and doing no measurement every other update of the 
whole space-time lattice, we have made of order 10,000,000 measurements, using blocks 
of 10,000 measurements to eliminate the effect of sticking 7 on the estimation of the 
statistical error. It takes almost one week of computing time to obtain one point in 
Fig.l with a processor executing 27 million floating-point operations per second. 8 

Our results and those of Moreo 4 suggest that it is the presence of strong 
antiferromagnetic correlations at half-filling which yield a maximum in the magnetic 
susceptibility away from half-filling. Fig. 3 of Ref. 6 shows that at U/t — 4 on 
a 4 x 4 lattice, the antiferromagnetic correlation length at half-filling has basically 
reached the size of the lattice 9 at the temperature we studied ((3 = 6/t), while it has 
not quite made it at the temperature {(5 = 4/t) studied by Moreo 4 . Evidently, at 
larger values of U/t, the antiferromagnetic correlation length could reach the system 
size even at f3 = 4/t. The results of Monte Carlo simulations on a related model also 
suggest the key role of strong antiferromagnetic correlations at half-filling: In the three- 



Figure 1: The static magnetic spin susceptibility Xs multiplied by temperature T is 
plotted as a function of band filling < n >. The units are t — 1, k B — 1, h — 2. Monte 
Carlo results are shown by points and error bars. The solid line is a guide to the eye. 
The dotted line shows the result for the free case (U=0) for the same system size and 
temperature. 

band model, Dopf et al. 10 have found a parameter range where the behavior observed 
experimentally is reproduced qualitatively and it seems that, in this parameter range, 
strong antiferromagnetic correlations are also present at half-filling. 11 

Whether these results of simulations on finite lattices have anything to do with 
the infinite system is a tricky question. A zeroth-order check consists in making sure 
that the U — result has the same qualitative dependence on filling for both the finite 
and the infinite lattice. In Fig.l, we plotted the U — result for a 4 x 4 lattice with a 
dotted line. Clearly, the spin susceptibility decreases as the system is doped away from 
half-filling, in qualitative agreement with the infinite system. This kind of agreement is 
not true at all fillings since on a finite lattice the spin susceptibility is not a monotonous 
function of filling while it is for the infinite system. Moreo 4 has argued that finite-size 



effects should be smaller for larger interaction strengths because of localization effects. 
This is almost certainly true, unless interactions introduce some other effect, such as 
a phase transition. We can guess what should happen at U/t = 4 by using some 
analytical results. It has been shown 12 that the magnetic structure factor obtained 
by Monte Carlo simulations is well described, not too close to half-filling, by the 
Generalized Random Phase Approximation (GRPA), if one takes into account two- 
particle correlations (Kanamori-Brueckner screening) by renormalizing the value of U. 
Using the renormalized value of U rn /t = 2.2 appropriate for the bare value U/t = 4, we 
find that long-range order sets in for a filling about equal to that where the maximum 
appears in Fig.l. This is consistent with the fact that in two- dimensions, the GRPA 
should start to fail when it predicts a magnetic phase transition. Indeed, a better 
approximation would take into account Mermin- Wagner fluctuations which prohibit 
long-range order in two- dimensions. How to include these fluctuations in itinerant 
electron theories is a problem which cannot presently be answered by Monte Carlo 
simulations but which is beginning to be successfully addressed analytically. 13 It would 
be reasonable to expect that in the regime where the mean-field phase transition 
found using the GRPA is suppressed by thermal fluctuations (T < T^ RPA ) the large 
antiferromagnetic fluctuations will nevertheless decrease the uniform spin susceptibility 
(at half-filling in Fig.l, it is smaller than the non-interacting value.). Because the 
tendency towards antiferromagnetism increases towards half-filling, the infinite-size 
two-dimensional Hubbard model would then show the same kind of behavior as the 
finite-size system of Fig.l. 14 Consistency with experiment might then occur not only 



in the strong-coupling limit, but also at relatively smaller couplings as long as Tq^ pa 
is not small compared with experimentally studied temperatures. In both the present 
work and that of Ref.4 however, the system sizes are much too small to yield a definitive 
answer to the question of the infinite-size limit. 
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